Prevalence and Mortality of Moose (Alces alces) Infected with Elaeophora schneideri in Montana, USA.

Elaeophora schneideri is a filarial nematode of North America that occasionally infects aberrant ruminant hosts such as moose (Alces alces). The role E. schneideri plays in clinical morbidity or mortality of moose remains uncertain. We sampled predominantly hunter-killed adult moose (n=127) to characterize the spatial patterns of prevalence and intensity of worms in carotid arteries of moose in Montana. We compared prevalence and intensity of E. schneideri within these moose to a separate sample of adult moose that died of health-related causes (n=34). We found lower prevalence in northwest Montana (0.06) than in the remainder of the state (0.42). We also found both higher prevalence of E. schneideri and higher intensity to be correlated with increased probability of health-related mortality. Our results suggest presence and intensity of E. schneideri correlate with mortality of moose, although the mechanisms of mortality remain uncertain.

Evaluating the summer landscapes of predation risk and forage quality for elk (Cervus canadensis).

The recovery of carnivore populations in North American has consequences for trophic interactions and population dynamics of prey. In addition to direct effects on prey populations through killing, predators can influence prey behavior by imposing the risk of predation. The mechanisms through which patterns of space use by predators are linked to behavioral response by prey and nonconsumptive effects on prey population dynamics are poorly understood. Our goal was to characterize population- and individual-level patterns of resource selection by elk (Cervus canadensis) in response to risk of wolves (Canis lupus) and mountain lions (Puma concolor) and evaluate potential nonconsumptive effects of these behavioral patterns. We tested the hypothesis that individual elk risk-avoidance behavior during summer would result in exposure to lower-quality forage and reduced body fat and pregnancy rates. First, we evaluated individuals' second-order and third-order resource selection with a used-available sampling design. At the population level, we found evidence for a positive relationship between second- and third-order selection and forage, and an interaction between forage quality and mountain lion risk such that the relative probability of use at low mountain lion risk increased with forage quality but decreased at high risk at both orders of selection. We found no evidence of a population-level trade-off between forage quality and wolf risk. However, we found substantial among-individual heterogeneity in resource selection patterns such that population-level patterns were potentially misleading. We found no evidence that the diversity of individual resource selection patterns varied predictably with available resources, or that patterns of individual risk-related resource selection translated into biologically meaningful changes in body fat or pregnancy rates. Our work highlights the importance of evaluating individual responses to predation risk and predator hunting technique when assessing responses to predators and suggests nonconsumptive effects are not operating at a population scale in this system.

Habitat selection by wolves and mountain lions during summer in western Montana.

In the Northern Rockies of the United States, predators like wolves (Canis lupus) and mountain lions (Puma concolor) have been implicated in fluctuations or declines in populations of game species like elk (Cervus canadensis) and mule deer (Odocoileus hemionus). In particular, local distributions of these predators may affect ungulate behavior, use of space, and dynamics. Our goal was to develop generalizable predictions of habitat selection by wolves and mountain lions across western Montana. We hypothesized both predator species would select habitat that maximized their chances of encountering and killing ungulates and that minimized their chances of encountering humans. We assessed habitat selection by these predators during summer using within-home range (3rd order) resource selection functions (RSFs) in multiple study areas throughout western Montana, and tested how generalizable RSF predictions were by applying them to out-of-sample telemetry data from separate study areas. Selection for vegetation cover-types varied substantially among wolves in different study areas. Nonetheless, our predictions of 3rd order selection by wolves were highly generalizable across different study areas. Wolves consistently selected simple topography where ungulate prey may be more susceptible to their cursorial hunting mode. Topographic features may serve as better proxies of predation risk by wolves than vegetation cover-types. Predictions of mountain lion distribution were less generalizable. Use of rugged terrain by mountain lions varied across ecosystem-types, likely because mountain lions targeted the habitats of different prey species in each study area. Our findings suggest that features that facilitate the hunting mode of a predator (i.e. simple topography for cursorial predators and hiding cover for stalking predators) may be more generalizable predictors of their habitat selection than features associated with local prey densities.

Multi-scale habitat assessment of pronghorn migration routes.

We studied the habitat selection of pronghorn (Antilocapra americana) during seasonal migration; an important period in an animal's annual cycle associated with broad-scale movements. We further decompose our understanding of migration habitat itself as the product of both broad- and fine-scale behavioral decisions and take a multi-scale approach to assess pronghorn spring and fall migration across the transboundary Northern Sagebrush Steppe region. We used a hierarchical habitat selection framework to assess a suite of natural and anthropogenic features that have been shown to influence selection patterns of pronghorn at both broad (migratory neighborhood) and fine (migratory pathway) scales. We then combined single-scale predictions into a scale-integrated step selection function (ISSF) map to assess its effectiveness in predicting migration route habitat. During spring, pronghorn selected for native grasslands, areas of high forage productivity (NDVI), and avoided human activity (i.e., roads and oil and natural gas wells). During fall, pronghorn selected for native grasslands, larger streams and rivers, and avoided roads. We detected avoidance of paved roads, unpaved roads, and wells at broad spatial scales, but no response to these features at fine scales. In other words, migratory pronghorn responded more strongly to anthropogenic features when selecting a broad neighborhood through which to migrate than when selecting individual steps along their migratory pathway. Our results demonstrate that scales of migratory route selection are hierarchically nested within each other from broader (second-order) to finer scales (third-order). In addition, we found other variables during particular migratory periods (i.e., native grasslands in spring) were selected for across scales indicating their importance for pronghorn. The mapping of ungulate migration habitat is a topic of high conservation relevance. In some applications, corridors are mapped according to telemetry location data from a sample of animals, with the assumption that the sample adequately represents habitat for the entire population. Our use of multi-scale modelling to predict resource selection during migration shows promise and may offer another relevant alternative for use in future conservation planning and land management decisions where telemetry-based sampling is unavailable or incomplete.

Functional responses in habitat selection: clarifying hypotheses and interpretations.

A fundamental challenge in habitat ecology and management is understanding the mechanisms generating animal distributions. Studies of habitat selection provide a lens into such mechanisms, but are often limited by unrealistic assumptions. For example, most studies assume that habitat selection is constant with respect to the availability of resources, such that habitat use remains proportional to availability. To the contrary, a growing body of work has shown the fallacy of this assumption, indicating that animals modify their behavior depending on the context at broader scales. This has been termed a functional response in habitat selection. Furthermore, a diversity of methods is employed to model functional responses in habitat selection, with little attention to how methodology might affect scientific and conservation conclusions. Here, we first review the conceptual and statistical foundations of methods currently used to model functional responses and clarify the ecological tests evaluated within each approach. We then use a combination of simulated and empirical data sets to evaluate the similarities and differences among approaches. Importantly, we identified multiple statistical issues with the most widely applied approaches to understand functional responses, including: (1) a complex and important role of random- or individual-level intercepts in adjusting individual-level regression coefficients as resource availability changes and (2) a sensitivity of results to poorly informed individual-level coefficients estimated for animals with low availability of a given resource. Consequently, we provide guidance on applying approaches that are insensitive to these issues with the goal of advancing our understanding of animal habitat ecology and management. Finally, we characterize the management implications of assuming similarity between the current approaches to model functional responses with two empirical data sets of federally threatened species: Canada lynx (Lynx canadensis) in the United States and woodland caribou (Rangifer tarandus caribou) in Canada. Collectively, our assessment helps clarify the similarities and differences among current approaches and, therefore, assists the integration of functional responses into the mainstream of habitat ecology and management.

Linking habitat selection and predation risk to spatial variation in survival.

A central assumption underlying the study of habitat selection is that selected habitats confer enhanced fitness. Unfortunately, this assumption is rarely tested, and in some systems, gradients of predation risk may more accurately characterize spatial variation in vital rates than gradients described by habitat selection studies. Here, we separately measured spatial patterns of both resource selection and predation risk and tested their relationships with a key demographic trait, adult female survival, for a threatened ungulate, woodland caribou (Rangifer tarandus caribou Gmelin). We also evaluated whether exposure to gradients in both predation risk and resource selection value was manifested temporally through instantaneous or seasonal effects on survival outcomes. We used Cox proportional hazards spatial survival modelling to assess the relative support for 5 selection- and risk-based definitions of habitat quality, as quantified by woodland caribou adult female survival. These hypotheses included scenarios in which selection ideally mirrored survival, risk entirely drove survival, non-ideal selection correlated with survival but with additive risk effects, an ecological trap with maladaptive selection and a non-spatial effect of annual variation in weather. Indeed, we found positive relationships between the predicted values of a resource selection function (RSF) and survival, yet subsequently incorporating an additional negative effect of predation risk greatly improved models further. This revealed a positive, but non-ideal relationship between selection and survival. Gradients in these covariates were also shown to affect individual survival probability at multiple temporal scales. Exposure to increased predation risk had a relatively instantaneous effect on survival outcomes, whereas variation in habitat suitability predicted by an RSF had both instantaneous and longer-term seasonal effects on survival. Predation risk was an additive source of hazard beyond that detected through selection alone, and woodland caribou selection thus was shown to be non-ideal. Furthermore, by combining spatial adult female survival models with herd-specific estimates of recruitment in matrix population models, we estimated a spatially explicit landscape of population growth predictions for this endangered species.

Preferred habitat and effective population size drive landscape genetic patterns in an endangered species.

Landscape genetics provides a framework for pinpointing environmental features that determine the important exchange of migrants among populations. These studies usually test the significance of environmental variables on gene flow, yet ignore one fundamental driver of genetic variation in small populations, effective population size, N(e). W(e) combined both approaches in evaluating genetic connectivity of a threatened ungulate, woodland caribou. We used least-cost paths to calculate matrices of resistance distance for landscape variables (preferred habitat, anthropogenic features and predation risk) and population-pairwise harmonic means of N(e), and correlated them with genetic distances, FST and D(c). Results showed that spatial configuration of preferred habitat and Ne were the two best predictors of genetic relationships. Additionally, controlling for the effect of Ne increased the strength of correlations of environmental variables with genetic distance, highlighting the significant underlying effect of Ne in modulating genetic drift and perceived spatial connectivity. We therefore have provided empirical support to emphasize preventing increased habitat loss and promoting population growth to ensure metapopulation viability.

Separating spatial search and efficiency rates as components of predation risk.

Predation risk is an important driver of ecosystems, and local spatial variation in risk can have population-level consequences by affecting multiple components of the predation process. I use resource selection and proportional hazard time-to-event modelling to assess the spatial drivers of two key components of risk--the search rate (i.e. aggregative response) and predation efficiency rate (i.e. functional response)--imposed by wolves (Canis lupus) in a multi-prey system. In my study area, both components of risk increased according to topographic variation, but anthropogenic features affected only the search rate. Predicted models of the cumulative hazard, or risk of a kill, underlying wolf search paths validated well with broad-scale variation in kill rates, suggesting that spatial hazard models provide a means of scaling up from local heterogeneity in predation risk to population-level dynamics in predator-prey systems. Additionally, I estimated an integrated model of relative spatial predation risk as the product of the search and efficiency rates, combining the distinct contributions of spatial heterogeneity to each component of risk.

Transcending scale dependence in identifying habitat with resource selection functions.

Multi-scale resource selection modeling is used to identify factors that limit species distributions across scales of space and time. This multi-scale nature of habitat suitability complicates the translation of inferences to single, spatial depictions of habitat required for conservation of species. We estimated resource selection functions (RSFs) across three scales for a threatened ungulate, woodland caribou (Rangifer tarandus caribou), with two objectives: (1) to infer the relative effects of two forms of anthropogenic disturbance (forestry and linear features) on woodland caribou distributions at multiple scales and (2) to estimate scale-integrated resource selection functions (SRSFs) that synthesize results across scales for management-oriented habitat suitability mapping. We found a previously undocumented scale-specific switch in woodland caribou response to two forms of anthropogenic disturbance. Caribou avoided forestry cut-blocks at broad scales according to first- and second-order RSFs and avoided linear features at fine scales according to third-order RSFs, corroborating predictions developed according to predator-mediated effects of each disturbance type. Additionally, a single SRSF validated as well as each of three single-scale RSFs when estimating habitat suitability across three different spatial scales of prediction. We demonstrate that a single SRSF can be applied to predict relative habitat suitability at both local and landscape scales in support of critical habitat identification and species recovery.

The role of translocation in recovery of woodland caribou populations.

Maintenance of viable populations of many endangered species will require conservation action in perpetuity. Efforts to conserve these species are more likely to be successful if their reliance on conservation actions is assessed at the population level. Woodland caribou (Rangifer tarandus caribou) were extirpated recently from Banff National Park, Canada, and translocations of caribou to Banff and neighboring Jasper National Park are being considered. We used population viability analysis to assess the relative need for and benefits from translocation of individuals among caribou populations. We measured stochastic growth rates and the probability of quasi extinction of four populations of woodland caribou with and without translocation. We used two vital rates in our analysis: mean adult female survival and mean number of calves per breeding-age female as estimates of mean fecundity. We isolated process variance for each vital rate. Our results suggested the Tonquin caribou population in Jasper is likely to remain viable without translocation, but that translocation is probably insufficient to prevent eventual extirpation of the two other populations in Jasper. Simulated reintroductions of caribou into Banff resulted in a 53-98% probability of >8 females remaining after 20 years, which suggests translocation may be an effective recovery tool for some caribou populations.

Resolving issues of imprecise and habitat-biased locations in ecological analyses using GPS telemetry data.

Global positioning system (GPS) technologies collect unprecedented volumes of animal location data, providing ever greater insight into animal behaviour. Despite a certain degree of inherent imprecision and bias in GPS locations, little synthesis regarding the predominant causes of these errors, their implications for ecological analysis or solutions exists. Terrestrial deployments report 37 per cent or less non-random data loss and location precision 30 m or less on average, with canopy closure having the predominant effect, and animal behaviour interacting with local habitat conditions to affect errors in unpredictable ways. Home-range estimates appear generally robust to contemporary levels of location imprecision and bias, whereas movement paths and inferences of habitat selection may readily become misleading. There is a critical need for greater understanding of the additive or compounding effects of location imprecision, fix-rate bias, and, in the case of resource selection, map error on ecological insights. Technological advances will help, but at present analysts have a suite of ad hoc statistical corrections and modelling approaches available-tools that vary greatly in analytical complexity and utility. The success of these solutions depends critically on understanding the error-inducing mechanisms, and the biggest gap in our current understanding involves species-specific behavioural effects on GPS performance.

Survival in the Rockies of an endangered hybrid swarm from diverged caribou (Rangifer tarandus) lineages.

In North America, caribou (Rangifer tarandus) experienced diversification in separate refugia before the last glacial maximum. Geographical isolation produced the barren-ground caribou (Rangifer tarandus groenlandicus) with its distinctive migratory habits, and the woodland caribou (Rangifer tarandus caribou), which has sedentary behaviour and is now in danger of extinction. Herein we report on the phylogenetics, population structure, and migratory habits of caribou in the Canadian Rockies, utilizing molecular and spatial data for 223 individuals. Mitochondrial DNA analyses show the occurrence of two highly diverged lineages; the Beringian-Eurasian and North American lineages, while microsatellite data reveal that present-day Rockies' caribou populations have resulted from interbreeding between these diverged lineages. An ice-free corridor at the end of the last glaciation likely allowed, for the first time, for barren-ground caribou to migrate from the North and overlap with woodland caribou expanding from the South. The lack of correlation between nuclear and mitochondrial data may indicate that different environmental forces, which might also include human-caused habitat loss and fragmentation, are currently reshaping the population structure of this postglacial hybrid swarm. Furthermore, spatial ecological data show evidence of pronounced migratory behaviour within the study area, and suggest that the probability of being migratory may be higher in individual caribou carrying a Beringian-Eurasian haplotype which is mainly associated with the barren-ground subspecies. Overall, our analyses reveal an intriguing example of postglacial mixing of diverged lineages. In a landscape that is changing due to climatic and human-mediated factors, an understanding of these dynamics, both past and present, is essential for management and conservation of these populations.